﻿UNIVERSITATEA “TRANSILVANIA” DIN BRAȘOV Laurențiu-Mihail IVANOVICI Procesarea Imaginilor îndrumar de laborator 2006 © 2003 EDITURA UNIVERSITĂȚII TRANSILVANIA BRAȘOV Adresa: 500030 Brasov, B-dul Eroilor, Nr 9 Tel/Fax: 0268 - 47 53 48 E-mail: editura@unitbv ro Tipărit la: Tipografia Universității “Transilvania” din Brașov B-dul Eroilor 9 tel/fax: 0268 47 53 48 Toate drepturile rezervate Editură acreditată de CNCSIS Adresa nr 1615 din 29 mai 2002 Referenți: Prof dr ing luliu Szekely Ș l dr ing Angel Cațaron Descrierea CIP a Bibliotecii Naționale a României IVANOVICI, LAURENȚIU MIHAIL Procesarea imaginilor: îndrumar de laborator / Laurențiu-Mihail Ivanovici - Brațov: Editura Universității “Transilvania”, 2006 Bibliogr ISBN (10) 973-635-674-4; ISBN (13) 978-973-635-674-2 004 932 Cuprins Cuprins i Mulțumiri v Cuvânt înainte ix 1 Noțiuni introductive de Qt 1 1 1 Clasa QApplication 1 1 2 Clasa Qlmage 1 1 2 1 Manipularea imaginilor 2 1 2 2 Atribute ale imaginii 2 1 2 3 Manipularea pixelilor 3 1 2 4 Formate de imagine 4 2 Noțiuni introductive 5 2 1 Percepția imaginilor 5 2 1 1 Structura ochiului uman Formarea imaginii 5 2 1 2 Lumina Luminanța Strălucirea 6 2 2 Modelul matematic al imaginii 7 2 2 1 Modelul continuu al imaginii 7 2 2 2 Modelul discret al imaginii 8 2 3 Ețantionarea imaginilor 8 2 4 Cuantizarea imaginilor 9 2 5 Imaginile digitale 9 3 îmbunătățirea imaginilor 15 3 1 Operațiile punctuale 15 3 2 Accentuarea contrastului 16 3 3 întinderea maxima a contrastului 17 3 4 Binarizarea imaginilor 18 3 5 Negativarea imaginilor 19 3 6 Decuparea imaginilor 20 4 Egalizarea histogramei 25 4 1 Histograma unei imagini 25 4 2 Egalizarea histogramei 26 4 2 1 Algoritmul de egalizare a histogramei 26 4 3 Observații 27 5 Transformari geometrice de bază 31 5 1 Translația 31 5 2 Rotația 31 5 3 Oglindirea 33 6 Zgomotul în imagini 37 6 1 Zgomotul cu distribuție uniforma 38 6 2 Zgomotul cu distribuție gaussiană 39 6 3 Zgomotul de tip “sare și piper" 41 6 4 Alte tipuri de zgomot 41 7 Filtrarea imaginilor 47 7 1 Filtrarea liniară a imaginilor 47 7 1 1 Filtrele de netezire 48 7 1 2 Filtrele trece-sus 51 7 2 Filtrarea neliniarăa a imaginilor 52 7 2 1 Filtrele de ordine 53 8 Transformari unitare 59 8 1 Transformari unitare bidimensionale 60 8 1 1 Proprietățile transformărilor unitare 62 8 2 Transformata Fourier discretăa 62 8 2 1 Transformata Fourier unidimensională 62 8 2 2 Transformarea Fourier bidimensională 63 8 3 Transformata cosinus discretaă 64 8 4 Transformata sinus discretăa 65 9 Compresia imaginilor 69 9 1 Clasificarea metodelor de compresie 69 9 2 Algoritmul Huffman 70 9 3 Algoritmul RLE 73 9 3 1 Algoritmul RLE pentru imagini binare 73 9 3 2 Algoritmul RLE pentru imagini în tonuri de gri 73 9 4 Compresia cu transformate 74 10 Segmentarea imaginilor 81 10 1 Segmentarea orientată pe regiuni 81 10 1 1 Prăaguirea histogramei 82 10 1 2 Segmentarea prin crețstere de regiuni 83 10 2 Segmentarea orientată pe contururi 83 10 2 1 Tehnicile de gradient 84 10 2 2 Operatorii compas 86 10 2 3 Identificarea trecerilor prin zero ale celei de-a doua derivate 87 Lista figurilor 90 Bibliografie 92 iii iv Dedicație Doresc sa dedic aceasta carte doamnei Delia Gârbacea v vi Mulțumiri "Țin să mulțumesc domnului profesor luliu Szekely, conducătorul catedrei de Electronică ți Calculatoare, din cadrul Universității “Transilvania” din Brașov, pentru încrederea pe care mi-a acordat-o prin încredințarea orelor de laborator de procesarea imaginilor Aț vrea să mulțumesc domnului profesor Vasile Buzuloiu, conducătorul Laboratorului de Analiza ți Prelucrarea Imaginilor (LAP), din cadrul Universității “POLITEHNICA” din Bucurețti, pentru sprijinul și îndrumarea acordate în formarea mea profesională Vreau să mulțumesc colegilor din cadrul LAPI, conferențiar Mihai Ciuc ți conferențtiar Constantin Vertan, pentru discuțtiile purtate de-a lungul anilor, în special cele care au vizat lucrăarile de laborator de procesarea imaginilor Nu în ultimul rand, vreau sa mulțumesc studenților care, cu sugestii sau comentarii, au contribuit la îmbunăatăațtirea conțtinutului acestui îndrumar de laborator vii Cuvânt înainte Procesarea imaginilor este un domeniu relativ recent, care evoluează rapid Aplicațiile sale se întâlnesc pretutindeni: în medicină, armată, industrie, artă sau acolo unde informația din mediul înconjurator este reprezentată sub forma de imagini Principala aplicație o reprezintă îmbunătațirea informației conțtinute de imagini în vederea interpretăarii de căatre un subiect uman sau pentru vederea artificială a roboților Odată cu răspandirea aparatelor de fotografiat digitale, procesarea imaginilor nu mai este domeniul exclusiv al oamenilor de țtiință sau al inginerilor Aplicațtiile de procesare a imaginilor au păatruns în casa oricăarui fotograf amator, fiind utilizate pentru retuțsarea fotografiilor digitale îndrumarul de față se adresează în special studenților de la facultățile tehnice, care urmează cursul de Procesarea imaginilor, dar ți celor care doresc săa îsți însuțseascăa algoritmii fundamentali de procesare a imaginilor Lucrările de laborator prezentate abordează operațiile punctuale, de vecina-tate țsi cele integrale, de la accentuarea contrastului panăa la transformaări unitare Ultimul capitol, referitor la segmentarea imaginilor, reprezintă o incursiune în domeniul analizei imaginilor Scopul lucrărilor de laborator este acela de a prezenta noțiunile fundamentale de procesare a imaginilor alb-negru țsi de a exemplifica folosind limbajul C++ modalități de implementare a algoritmilor de procesare prezentați Pentru o înțtelegere cat mai bunăa a soluțtiilor propuse, cititorul trebuie săa fie familiarizat cu limbajele de programare C ți C++ ix Lucrarea 1 Noțiuni introductive de Qt Qt1 este un mediu de dezvoltare care include biblioteci de clase C++ și unelte de dezvoltare de aplicații pentru diverse sisteme de operare: Microsoft Windows, MAC OS sau Linux Bibliotecile Qt cuprind peste 400 de clase C++, care încapsulează infrastructura necesară pentru dezvoltarea de aplicații Qt API1 2 include clase pentru realizarea de interfețe grafice utilizator, pentru programarea în rețea, baze de date sau integrare OpenGL3 1 1 Clasa QApplication Clasa QApplication controleaza ți gestionează interfața grafică a unei aplica-țtii Conțtine bucla principalăa de evenimente a aplicațtiei, în care sunt procesate evenimentele provenind de la sistemul de ferestre De asemenea ges-tioneazăa configurăarile aplicațtiei, fazele de inițtializare țsi terminare a aplicațtiei Orice aplicație Qt este un obiect QApplication indiferent de numarul de ferestre ale aplicației Obiectul este accesibil prin apelarea funcției membre instance(), care returneaza un pointer către acest obiect 1 2 Clasa Qlmage Clasa Qlmage pune la dispoziție o reprezentare independentă de platformă a unei imagini, care permite accesul direct la valorile pixelilor Qt oferaă patru clase pentru manipularea imaginilor: Qlmage, QPixmap, QBitmap ți QPicture Clasa Qlmage este proiectata ți optimizata pentru operații de intrare/ieți-re ți pentru acces direct la pixelii imaginii, pe cand clasa QPixmap este proiectată ți optimizata pentru vizualizarea imaginilor pe ecran Clasa QBitmap 1 Trolltech, http://www trolltech com 2Application Program Interface 3 Limbaj de grafică 3D dezvoltat de firma Silicon Graphics 1 2 LUCRAREA 1 NOȚIUNI INTRODUCTIVE DE QT moștenește clasa Pixmap și este folosită pentru imagini binare, iar clasa QPicture este un “paint device" în continuare ne vom concentra atenȘia asupra clasei QImage Clasa QImage poate manipula o serie de formate de imagine, care includ fie imagini monocrome, fie reprezentate pe 8 biȘi sau 32 biȘi Clasa pune la dispoziȘie o colecȘie de funcȘii care pot fi folosite pentru obȘinerea de informaȘii despre imagine sau care permit anumite transformări ale imaginii 1 2 1 Manipularea imaginilor Clasa QImage ofera cateva moduri de a încarca o imagine dintr-un fiȘier: la instanȘierea obiectului Qîmage sau folosind funcȘiile loadFromData() sau load(), care vor fi apelate după crearea obiectului Qlmage Pentru a salva un obiect de tip Qlmage se foloseȘte funcȘia save() Următorul exemplu ilustreaza folosirea funcȘiei load(): Qlmage image; QString filename = QFileDialog :: getOpenFileName(); image load( filename, 0 ); Lista completa a formatelor de fiȘiere cunoscute este disponibilă prin intermediul a două metode: QImageReader::supportedImageFormats() Și QImageWriter::supportedImageFormats() Implicit, Qt suportă urmatoa-rele formate de fiȘsiere: Format Descriere Operații BMP Windows Bitmap Read/Write GIF Graphic înterchange Format (opȘional) Read JPG Joint Photographic Experts Group Read/Write JPEG Joint Photographic Experts Group Read/Write PNG Portable Network Graphics Read/Write PBM Portable Bitmap Read PGM Portable Graymap Read PPM Portable Pixmap Read/Write XBM X11 Bitmap Read/Write XPM X11 Pixmap Read/Write Tabelul 1 1: Formate de fișiere imagine suportate de Qt 1 2 2 Atribute ale imaginii QImage oferăa o colecȘtie de funcȘtii pentru obȘtinerea de informaȘtii despre atributele imaginii, cum ar fi geometria imaginii sau informaȘtii despre paleta de culori 1 2 CLASA QIMAGE 3 Atribute Funcții Geometrie Funcțiile size(), width(), height(), dotsPerMeterX() ți dotsPerMeterY() oferaă informațtii despre dimensiunile țsi aspectul imaginii Culoare Culoarea unui pixel poate fi aflatăa spe-cificand coordonatele pixelului funcției pixel(), care returneazaă culoarea ca o valoare QRbg In cazul imagini- lor monocrome sau cu nivele de gri, funcțtiile numColors() țsi colorTable oferăa informațtii despre componentele de culoare ale imaginii Nivel jos Funcția depth() returnează numărul de bițti pe care este reprezentataă valoarea unui pixel: 1 (imagini monocrome), 8 sau 32 Funcțiile format(), bytesPerLine() sți numBytes() oferăa informațtii despre modul de stocare a imaginii Tabelul 1 2: Funcții de aflare a atributelor unei imagini 1 2 3 Manipularea pixelilor în cazul în care culorile pixelilor sunt stocate pe 32 biți, funcția setPixel se poate utiliza pentru modificarea valorii unui pixel specificat prin coordonatele sale Noua sa valoare este un cvadruplu ARGB4, fiind unul din argumentele funcției Pentru a specifica valoarea unui pixel în modelul RGB se folosețte funcția qRgb care returnează un obiect QRgb In continuare este prezentat un exemplu de folosire a funcției setPixel Se crează mai întâi un obiect QImage, reprezentand o imagine de dimensiune 3x3 QImage image(3, 3, QImage::Format RGB32); QRgb value; value = qRgb(189, 149, 39); // 0xffbd9527 image setPixel(1, 1, value); value = qRgb(122, 163, 39); // 0xff7aa327 image setPixel(0, 1, value); image setPixel(1, 0, value); 4Alpha Red Green Blue 4 LUCRAREA 1 NOȚIUNI INTRODUCTIVE DE QT value = qRgb(237, 187, 51); // 0xffedba31 image setPixel(2, 1, value); 1 2 4 Formate de imagine Fiecare pixel stocat într-un obiect de tip QImage este reprezentat ca un întreg Dimensiunea acestui întreg depinde de format Qt suporta imagini monocrome, reprezentate pe 1 bit și imagini reprezentate pe 8 sau 32 biți Imaginile monocrome sunt stocate folosind indici de 1 bit într-o tabelă de culoare cu doar doua intrări (culori) In funcție de ordinea de stocare a biților, big endian sau little endian, diferențiem două tipuri de imagini monocrome Imaginile reprezentate pe 8 bițti sunt stocate folosind indici de 8 bițti într-o tabelă de culoare, avand un byte per pixel Tabela de culoare este un obiect QVector Imaginile reprezentate pe 32 biți nu au tabela de culoare, fiecare pixel conținand o valoare QRgb Cea mai râspandită modalitate de a reprezenta tripletul RGB pe 32 biți este urmatoarea: GxffRRGGBB, în care se folosesc cate 8 biți pentru fiecare componentă Formatul unei imagini se poate determina folosind funcția format() Pentru a converti o imagine într-un nou format se folosețte funcția convertToFormat(), care acceptăa ca argument noul format al imaginii Lucrarea 2 Noțiuni introductive de prelucrarea imaginilor BREVIAR TEORETIC 2 1 Percepția imaginilor înțelegerea procesului de percepție a imaginilor de către ochiul uman este foarte importantă pentru dezvoltarea de tehnici de evaluare a calitații unui sistem sau algoritm de procesare a imaginilor Informația la nivel vizual, conținută de către o imagine, reprezintă o distribuție spațială a unei anumite mărimi, cum ar fi de exemplu luminanța obiectelor ce compun respectiva imagine Informația percepută de ochiul uman poate fi definita de atribute cum ar fi strălucirea, culoarea sau muchiile unui obiect 2 1 1 Structura ochiului uman Formarea imaginii Ochiul uman are o formă aproape sferică, avand un diametru, în medie, de aproximativ 2 cm Este format din mai multe membrane: corneea si sclerotica, ca învelis exterior, ți coroida ți retina în interior Corneea este un țesut dur si transparent, pe cand sclerotica este o țesut opac Coroida se află imediat sub sclerotica ți conține o rețea de vase de sange, ce hrănesc ochiul Este puternic pigmentataă pentru a reduce excedentul de luminaă ce intrăa în ochi Pe cea mai din interior membrană, retina, se formează imaginea, sub influențta luminii reflectate de obiectele exterioare ochiului Retina conține două tipuri de receptori: conuri și bastonade Conurile, în număr de aproximativ 6-7 milioane, servesc la percepția culorilor Vederea umana poate percepe detalii foarte fine datorită densității mari a acestor receptori, fiecare con fiind legat la o terminație nervoasă Bastonațele, în număr de aproximativ 75-150 milioane, servesc vederii crepusculare, în condiții de iluminare slabă Acestea sunt răspandite pe o arie mare ți conec 5 6 LUCRAREA 2 NOȚIUNI INTRODUCTIVE tate la doar câteva terminații nervoase, ceea ce are ca efect reducerea considerabilă a perceperii detaliilor din imagine Cristalinul joaca rolul de lentila El conține aproximativ 70% apa ți, prin compoziția lui, permite trecerea a doar 8% din spectrul de radiație vizibilă, absorbind în bună măsura ți radiațiile infraroții ți ultraviolete Imaginea obiectelor se proiectează pe retină, prin cristalin, fiind râsturna-tăa țsi avand dimensiuni mult mai mici, principiu ce staă la baza aparatului de fotografiat (vezi Figura 2 1) Figura 2 1: Formarea imaginii în aparatul de fotografiat: 1-lentilă, 2-camera obscură, 3-obiectiv, 4-peliculă 2 1 2 Lumina Luminanța Strălucirea Lumina este o radiație electromagnetică ce stimulează receptorii de la nivelul retinei Ea se exprimă ca fiind o distribuție L(A) de energie, unde A este lungimea de undă a radiației, în cazul nostru vizibilă, cu valori între 350 ți 780 nm Lumina percepută de la un obiect se poate scrie matematic astfel: I (A) = p(A)L(A) (2 1) unde p(A) reprezinta măsura în care un obiect reflectă sau transmite energia luminoasa incidentă, a carei distribuție este exprimată prin L(A) Luminanța sau intensitatea luminoasa a unui obiect cu o distribuție spațială a luminii, I(x,y,A), se definețte astfel: f (x, У) = f I(x, У, A)V (A)dA (2 2) 0 unde V(A) este funcția de eficiență luminoasă relativă a sistemului vizual Pentru ochiul uman, V(A) este o curbă de tip clopot, a cărei caracteristici depind de la o persoanăa la alta (vezi Figura 2 2) 2 2 MODELUL MATEMATIC AL IMAGINII 7 Funcția de eficienta luminoasa relativa 400 0 440 0 480 0 520 0 560 0 600 0 640 0 680 0 720 0 760 0 —I—’—I—t—I—’—I—t—I—’—I—t—I—t—I—t—I—’—I— Eficienta relativa 400 0 440 0 480 0 520 0 560 0 600 0 640 0 680 0 720 0 760 0 lambda (nm) Figura 2 2: Forma tipică a funcției de eficiență luminoasă relativă Luminanța unui obiect este independentă de luminanța obiectelor din jur Strălucirea unui obiect este luminanța percepută și depinde de luminanța mediului ambiant obiectului Două obiecte aflate într-un același ambient pot avea aceeași luminanța, dar străluciri diferite 2 2 Modelul matematic al imaginii 2 2 1 Modelul continuu al imaginii Matematic imaginile pot fi reprezentate ca o funcție de două variabile, în spațiul L2(E2), astfel: • imaginile în tonuri de gri se pot modela cu: f(x,y) : R2 R, caz în care valorile funcției / reprezintă valorile luminanței obiectelor din imagine, în punctele (x, y) ale spațiului • imaginile color se pot modela cu: f(x,y) : R2 —>■ R3, caz în care valorile funcției / reprezintă vectori de 3 componente dintr-un spațiu al culorilor De exemplu pot fi cele trei componente ale modelului RGB1 Spațiul L2(R2) poate fi limitat la un domeniu finit D, ca în Figura 2 3 XRGB = (engl ) Red Green Blue 8 LUCRAREA 2 NOȚIUNI INTRODUCTIVE Figura 2 3: Domeniu finit din R2 2 2 2 Modelul discret al imaginii Acesta este modelul utilizat în practică Funcția f ia valori discrete, fiind deasemenea definită pe un domeniu de valori discrete, adică: f (k, l) : Z2 ! Z+ sau f (k, l) : Z2 ! Z+3 (2 3) Trecerea de la domeniul continuu la domeniul discret se face prin eșantionare ți cuantizare 2 3 Eșantionarea imaginilor Pentru a putea prelucra cu ajutorul unui calculator o imagine f (x, y), aceasta trebuie discretizată spațial și în amplitudine Discretizarea coordonatelor spațiale (x, y) poartă numele de eșantionare Eșantionarea reprezintă procesul de aproximare a unei imagini continue f (x,y) cu o matrice de dimensiune MxN, astfel: f (x,y) f (0, 0) f(1, 0) f (0,1) f (0,M - 1) \ f(1,1) f(1,M - 1) (2 4) \f (N - 1, 0) f (N - 1,1) f (N - 1,M - 1) / Teorema eșantionării: O imagine f(x,y), avand un spectru finit, ețantionata uniform cu o rețea dreptunghiulară de forma celei din Figura 2 4 poate fi refăcuta fără eroare din ețantioanele f (mâx,ndy) cu ajutorul formulei de interpolare: sin^(xus — m) sin^(xvs — n Ux», — m) %(xvs — n) (2 5) 2 4 CUANTIZAREA IMAGINILOR 9 unde us și vs reprezintă frecvențele spațiale de eșantionare Egalitatea dată de teorema eșantionării este valabilă dacă ți numai dacă este respectată condiția Nyquist, ți anume: 1 — = us > 2uo, ox 1 — = vs > 2vo Oy (2 6) unde u0 ți v0 reprezinta frecvențele spațiale maxime care apar în imagine n dx dy m Figura 2 4: Rețea dreptunghiulară de ețantionare Cu alte cuvinte, frecvențele spațiale de ețantionare trebuie să fie cel puțin dublul frecvențtelor spațtiale maxime conțtinute de imagine 2 4 Cuantizarea imaginilor Cuantizarea este procesul de discretizare a valorilor funcției f (x, y) Aceasta se realizează de obicei cu ajutorul unei funcții de tip scară, de forma celei din Figura 2 5 Astfel, tuturor valorilor lui x dintr-un interval li se atribuite valori discrete k Cuantizarea este un proces însoțit de zgomot, cunoscut sub numele de eroare de cuantizare Cea mai utilizata metodă de cuantizare este cea uniforma, ceea ce înseamna că intervalele funcției de cuantizare sunt egale 2 5 Imaginile digitale Imaginile astfel discretizate reprezintă structuri bidimensionale de date, denumite imagini digitale Un element (k,l) al imaginii poartă numele de pixel2 2pixel = (engl picture + element) 10 LUCRAREA 2 NOȚIUNI INTRODUCTIVE k Figura 2 5: Exemplu de funție de cuantizare Imaginile digitale pot fi stocate în memoria sau pe discul unui sistem de procesare ți analiză a imaginilor, în vederea vizualizării sau prelucrării ulterioare Pe disc imaginile sunt stocate sub forma unor fișiere Fișierele pot fi de mai multe feluri, în funcție de formatul în care pastreaza datele ce reprezinta imagini: BMP, JPEG, GIF, TIFF, etc Cel mai simplu format de fișier imagine este Windows Bitmap (BMP) al firmei Microsoft Acesta este formatul în care o imagine digitală este stocată practic fără nici un fel de codare sau pierdere de informație, cu excepția reprezentarii binare DESFĂȘURAREA LUCRĂRII In continuare este prezentat codul C++ al unei aplicați care citețte o imagine în format BMP ți o afisează într-o fereastră pe ecran Aplicația este dezvoltată utilizand bibliotecile Qt (Linux) Citiți ți înțelegeți codul Fițierul aplicatie h conține declararea clasei ImageViewer #ifndef APLICATIE H #define APLICATIE H #include #include #include #include class QMenuBar; class QPopupMenu; class ImageViewer : public QWidget 2 5 IMAGINILE DIGITALE 11 { Q OBJECT public: ImageViewer( QWidget *parent = 0, const char *name = 0, int wFlags = 0 ); ~ImageViewer(); bool loadImage( const char *fileName ); protected: void paintEvent( QPaintEvent * ); private: int conversion flags; int alloc context; QImage image; QPixmap pm; QMenuBar *menubar; QPopupMenu *file; QLabel *status; bool reconvertImage(); private slots: void openFile( void ); void saveFile( void ); }; #endif // APLICATIE H Fișierul aplicatie cpp conține definițiile funcțiilor membre clasei Im-ageViewer #include "aplicatie h" #include #include #include #include #include #include #include ImageViewer :: ImageViewer( QWidget *parent, const char *name, int wFlags ) : QWidget( parent , name , wFlags ), conversion flags( PreferDither ), filename( 0 ) { alloc context = 0; 12 LUCRAREA 2 NOȚIUNI INTRODUCTIVE menubar = new QMenuBar( this ); menubar->setSeparator( QMenuBar :: InWindowsStyle ); file = new QPopupMenu(); CHECK PTR( file ); menubar->insertItem( "&File" , file ); file->insertItem( "Open image ( BMP )", this, SLOT( openFile() ), CTRL+Key O ); file->insertSeparator(); file->insertItem( "Save as ( BMP )", this, SLOT( saveFile() ), CTRL+Key S ); file->insertSeparator(); file->insertItem( "Quit", qApp, SLOT(quit()), CTRL+Key Q ); status = new QLabel( this ); status->setFrameStyle( QFrame::WinPanel | QFrame::Sunken ); status->setFixedHeight( fontMetrics() height() + 4 ); setMouseTracking( TRUE ); } ImageViewer :: ~ImageViewer() { if ( alloc context ) QColor :: destroyAllocContext( alloc context ); delete menubar; delete file; delete help; delete status; } void ImageViewer :: openFile() { QString newfilename = QFileDialog :: getOpenFileName(); if ( !newfilename isEmpty() ) { loadImage( newfilename ); repaint(); } } void ImageViewer :: saveFile( void ) 2 5 IMAGINILE DIGITALE 13 { QImageIO iio; QString save as file = QFileDialog :: getSaveFileName(); if( !save as file isEmpty() ) { iio setImage( image ); iio setFileName( save as file ); iio setFormat( "BMP" ); iio write(); } } bool ImageViewer :: loadImage( const char *fileName ) { bool ok = FALSE; int w, h; if( filename ) { QApplication :: setOverrideCursor( waitCursor ); ok = image load( filename , 0 ); pm convertFromImage( image , conversion flags ); if( ok ) { setCaption( filename ); w = pm width(); h = pm height(); h += menubar->heightForWidth(w) + status->height(); } else { pm resize( 0,0); update(); } setFixedWidth( w ); setFixedHeight( h ); status->setGeometry( 0, height() - status->height(), width(), status->height() ); QApplication :: restoreOverrideCursor(); } return ok; 14 LUCRAREA 2 NOȚIUNI INTRODUCTIVE } void ImageViewer :: paintEvent( QPaintEvent *e ) { if( pm size() != QSize( 0,0)) { QPainter painter( this ); painter setClipRect( e->rect() ); painter drawPixmap( 0, menubar->heightForWidth( width() ), pm ); } } Fișierul main cpp conține instanțierea unui obiect QApplication și a unuia de tip ImageViewer #include "aplicatie h" #include #include #ifdef QIMGIO #include #endif int main( int argc, char **argv ) { QApplication :: setFont( QFont( "Helvetica" , 12 ) ); QApplication a( argc, argv ); #ifdef QIMGIO qInitImageIO(); #endif ImageViewer *w = new ImageViewer( 0, "main window", QWidget :: WDestructiveClose ) ; w->show() ; QObject :: connect( qApp, SIGNAL( lastWindowClosed() ), qApp, SLOT(quit()) ); return a exec() ; } Lucrarea 3 îmbunătățirea imaginilor prin operații punctuale BREVIAR TEORETIC Termenul general de îmbunătățire a imaginilor se referă la o clasă largă de operații, ce au ca scop mărirea detectabilității componentelor imaginii Această detectabilitate depinde în primul rând de percepția vizuală a unui observator uman ți deci reprezinta o apreciere subiectivă a imaginii îmbunătațirea calității unei imagini se face făra a presupune vreun model de degradare sau a lua în considerare vreo informațtie legatăa de imaginea originală Paradoxal, dar ți o imagine originală (nedegradată) poate fi îmbunaătăațtitaă în general, îmbunaătăațtirea se referăa la accentuarea unor caracteristici ale imaginii, cum ar fi muchiile, contururile sau contrastul Procesul de îmbunaătăațtire nu măarețste cantitatea de informațtie conțtinutăa de o imagine 3 1 Operațiile punctuale Operațiile punctuale sunt definite de o funcție, care atribuie un nou nivel de gri pixelilor din imagine Noua valoare a pixelului va depinde doar de vechea valoare a acestuia, de unde ți denumirea de “operație punctuală" Matematic se poate scrie: g(k,l) = ^(/(k,l)) (3 1) unde f este imaginea originală, iar g imaginea îmbunatățita g(k,l) reprezintă noua valoare a pixelului (k, l), iar f (k, l) vechea valoare Operația punctuală este descrisă de funcția ф O operație sau o transformare punctuală este reprezentată în Figura 3 1 15 16 LUCRAREA 3 ÎMBUNĂTĂȚIREA IMAGINILOR Figura 3 1: Reprezentarea unei operații punctuale 3 2 Accentuarea contrastului Această operație se folosește în special pentru îmbunătațirea imaginilor cu contrast scăzut Acesta poate apărea datorită unei slabe iluminări, a unei iluminări neuniforme sau din cauza unor neliniarități ale caracteristicilor unui senzor de captură a imaginii O funcție tipică de accentuare a contrastului are urmăatoarea formaă matematicăa: = T unde T este o valoare de prag, reprezentand o valoare întreagă din intervalul [0; L) Pentru imaginile color, binarizarea se poate face în urmăatorul mod: Ф(Ѵ) = 0; L - 1; pentru Y(v) T (3 5) unde v este un vector tridimensional ce reprezintaă culoarea pixelului (de exemplu v = (R;G;B) ), iar Y(v) reprezintă luminanța (Y = 0 3R + 0 6G + 0 1B ) In urma acestei transformări, contrastul este maximizat la nivelul întregii imagini 3 5 NEGATIVAREA IMAGINILOR 19 0(хф L-1 0 T L-1 x Figura 3 5: Funcția de binarizare (a) (b) Figura 3 6: Binarizarea: (a) imaginea originală și (b) imaginea binarizata cu T=127 3 5 Negativarea imaginilor Negativul unei imagini se obține prin inversarea ordinii nivelelor de gri Pentru imaginile în tonuri de gri, operația de negativare se face folosind funcția: ф(х) = (L — 1) — x reprezentata în Figura 3 7, iar pentru imaginile color: ф(ѵ) = ((L - 1) - R; (L - 1) - G; (L - 1) - B) (3 6) (3 7) 20 LUCRAREA 3 ÎMBUNĂTĂȚIREA IMAGINILOR Figura 3 7: Funcția de negativare Negativarea imaginilor este utilă pentru afișarea unor imagini medicale sau pentru realizarea de imagini negative pe suporturi fizice, de exemplu, de tip peliculăa (a) (b) Figura 3 8: Negativarea: (a) imaginea originala și (b) negativul imaginii 3 6 Decuparea imaginilor Decuparea cu păstrarea fundalului este data de formula: ^(x) = L - 1 x, pentru x 2 [a,b] în rest (3 8) iar decuparea faărăa păastrarea fundalului este datăa de formula: 3 6 DECUPAREA IMAGINILOR 21 L - 1; 0; pentru x 2 [a;b] în rest (3 9) Figura 3 9: Funcția de decupare cu păstrarea fundalului 0(X) L-1 0 a b L-1 x Figura 3 10: Funcția de decupare fără păstrarea fundalului Aceste operații permit “decuparea” unor regiuni din imagine, caracterizate de anumite nivele de gri Acest lucru este folosit atunci cand diferite carcaterisitici ale imaginii sînt conținute în nivelele de gri respective, cum ar fi de exemplu decuparea regiunilor de temperatură joasa reprezentate de nori din imaginile obținute de un satelit meteo în astfel de imagini, nivelele de gri ce corespund unor nori sînt direct proporțtionale cu valorile de temperaturi joase DESFĂȘURAREA LUCRĂRII Problema 1 Compilați sursele C++ ale lucrarii Rulați aplicația ți 22 LUCRAREA 3 ÎMBUNĂTĂȚIREA IMAGINILOR observați rezultatul accentuării contrastului pentru o imagine în tonuri de gri (lena AN bmp) Funcția care realizează accentuarea contrastului este prezentată în continuare, pentru urmatoarele valori: a = 80, b = 170, Va = 20 ți Vb = 235 Pantele care caracterizează funcția de accentuare a contrastului vor fi urmatoarele: — 255-170 — 85 ■ int ImageViewer :: f accentuare( int nivel gri ) { if( nivel gri >= 0 && nivel gri 80 && nivel gri 170 && nivel gri = w ) iprim -= w; 5 3 OGLINDIREA 35 imag tx setPixel( iprim, j, pixel ); } image = imag tx; pm = image; update(); } Problema 1 Implementați operația de translație, pentru următorii doi vectori de translație: (0,100) ți (100,100) Problema 2 Implementați operația de oglindire “stânga-dreapta” față de axa verticală ce trece prin centrul imaginii Problema 3 Implementați operația de oglindire “sus-jos” față de axa orizontală ce trece prin centrul imaginii 36 LUCRAREA 5 TRANSFORMĂRI GEOMETRICE DE BAZA Lucrarea 6 Zgomotul în imagini BREVIAR TEORETIC Zgomotul este un semnal aleator, care afectează informația utilă conținută într-o imagine El poate apare de-alungul unui lanț de transmisiune, sau prin codarea ți decodarea imaginii, ți reprezintă un element perturbator nedorit De obicei se încearca eliminarea lui prin diverse metode de filtrare Zgomotul se poate suprapune informației utile în două moduri: • aditiv în acest caz, zgomotul se numețte zgomot aditiv ți matematic se poate scrie: g(x,y) = f (x,y) + n(x,y) (6 1) unde f (x,y) este imaginea inițiala, neafectata de zgomot, n(x,y) este zgomotul, iar g(x, y) este imaginea afectată de zgomot • multiplicativ în acest caz zgomotul se numețte zgomot multiplicativ, iar fenomenul de suprapunere al acestuia peste informațtia utilaă se scrie matematic astfel: g(x,y)= f (x,y) ■ n(x,y) (6 2) unde f (x, y), n(x, y) ți g(x, y) au aceleați semnificații ca mai sus în continuare vom trata zgomotul ca fiind aditiv Modelul de degradare a imaginii este reprezentat în figura 6 1 Zgomotul multiplicativ poate fi tratat la fel de simplu ca zgomotul aditiv, dacăa logaritmaăm relațtia (6 2): log(g(x, y)) = log(f (x, y) ■ n(x, y)) = log(f (x, y)) + log(n(x, y)) (6 3) Cantitativ, se poate aprecia măsura în care zgomotul a afectat informația utila, calculand raportul semnal-zgomot1: 1SNR = (engl ) Signal Noise Ratio 37 38 LUCRAREA 6 ZGOMOTUL ÎN IMAGINI f(x,y) n(x,y) g(x,y) Figura 6 1: Modelul aditiv de degradare M-1 N-1 E Ef2 (i;j) SNR = 10log -[dB ] (6 4) Ё Ё n2(i,j) i=0 j=0 M-1 N-1 E E f2(j SNR = 10logM 1 N=0 j=0 -[dB] (6 5) Ё Ё [f (j - g(i;j)]2 i=0 j=0 Raportul semnal-zgomot reprezintă raportul dintre energia imaginii originale și energia zgomotului suprapus acesteia în continuare vom asocia valorilor pe care le ia zgomotul o variabilă aleatoare £, care, în funcție de tipul zgomotului, va fi caracterizată de diferite funcții de densitate de probabilitate 6 1 Zgomotul cu distribuție uniformă Zgomotul cu distribuție uniformă (vezi Figura 6 2) este caracterizat de o funcțtie de densitate de probabilitate de forma: (x) = A, 0; pentru x 2 [—2; 2] , în rest (6 6) In Figura 6 3 puteți observa efectele zgomotului cu distribuție uniformă asupra imaginii “Lena”, pentru un raport semnal/zgomot de 5 dB 6 2 ZGOMOTUL CU DISTRIBUȚIE GAUSSIANA 39 k/2 w(x) A 0 k/2 x Figura 6 2: Funcția de densitate de probabilitate pentru o distribuție uniformă (a) (b) Figura 6 3: Zgomotul uniform: (a) imaginea originală; (b) imaginea afectată de zgomot uniform, SNR=5 dB 6 2 Zgomotul cu distribuție gaussiană Zgomotul gaussian este caracterizat de o funcțtie de densitate de probabilitate de forma (vezi Figura 6 4): 1 (x-M)2 wț(x) = V 2%a2 (6 7) e 2^'2 40 LUCRAREA 6 ZGOMOTUL IN IMAGINI Figura 6 4: Funcția de densitate de probabilitate pentru o distribuție gaussiană Valoarea medie a unei variabilei aleatoare £, cu o distribuție dată de funcția wț(x) ca în relația (6 7), este /л, iar varianța ei este ст2 O distribuție gaussiană se mai numește și normală O distribuție normală de medie p, și varianță ст2 se notează cu V(/u,ct2) în Figura 6 5 puteți observa efectele zgomotului cu distribuție gaussiană asupra imaginii “Lena”, pentru un raport semnal/zgomot de 5 dB (a) (b) Figura 6 5: Zgomotul Gaussian: (a) imaginea originală; (b) imaginea afectată de zgomot gaussian, SNR=5 dB 6 3 ZGOMOTUL DE TIP “SARE ȘI PIPER" 41 6 3 Zgomotul de tip “sare si piper" După cum îi spune numele, acest tip de zgomot va afecta valorile pixelilor în doua moduri: “sare” - adica noua valoare a pixelului va fi 255 (pixelul va fi alb), sau “piper” - adică noua valoare a pixelului va fi 0 (pixelul va fi negru) Zgomotul de tip “sare și piper” (vezi Figura 6 6) este perfect decorelat, deoarece între valorile 0 și 255, și între coordonatele pixelilor afectați de zgomot nu exista corelaȘie (a) (b) Figura 6 6: Zgomotul “sare și piper”: (a) imaginea originala; (b) imaginea cu 10% pixeli afectați de zgomot 6 4 Alte tipuri de zgomot Zgomotele diferaă între ele în funcștie de distribuștia care le caracterizeazăa Alte funcții de distribuție utilizate sunt: • distribuția Rayleigh: (x) = xe 2 • distribuția Maxwell: w^(x) = x2e • distribuția Beta: w^ (x) = xb(1 — x)c • distribuția Gamma (Erlang): w^ (x) = xne~x • distribuția Laplace: w^(x) = e jxj • distribuția Cauchy: w^ (x) = -+2 Aceste funcții sunt reprezentate în Figura 6 7 42 LUCRAREA 6 ZGOMOTUL IN IMAGINI Figura 6 7: Diverse funcții de distribuție DESFĂȘURAREA LUCRĂRII Problema 1 Suprapuneți zgomot cu distribuție uniformă peste o imagine în tonuri de gri (în meniul Algoritmi, funcția zgomot uniform) 6 4 ALTE TIPURI DE ZGOMOT 43 void ImageViewer :: zgomot uniform( void ) { int i, j; int w = image width(); int h = image height(); long int e zgomot = 0; //energia zgomotului long int e imagine = 0; //energia imaginii double SNR; //raportul semnal-zgomot //imaginea afectata de zgomot QImage image aff( w, h, 32, 0, QImage::IgnoreEndian ); int med = 0; //media zgomotului int dis = 200; //dispersia zgomotului for( i = 0; i 255 ) val = 255; if( val = 50 ) sare piper = qRgb( 255, 255, 255 ); else sare piper = qRgb( 0, 0, 0 ); if( (i >= 0) && (i = 0) && (j 0 Relațtia (7 4) poartaă numele de condițtie de normare pentru filtre de netezire (sau trece-jos) Filtrul de mediere Filtrul de mediere este cel mai simplu filtru de netezire Caracteristic unui filtru de mediere este faptul ca toți coeficienții măștii de filtare sânt egali Dacă ținem cont ți de condiția de normare, atunci coeficienții unui filtru de mediere, care are o masca de filtrare de dimensiune N x N, au valoarea N2 ân Figurile 7 2, 7 3 ți 7 4 sunt prezentate mățtile de filtrare de mediere, pentru N = 3, 5 ți respectiv 7 1 1 1 9 9 9 1 1 9 1 1 9 9 Figura 7 2: Masca de filtrare patrată de dimensiune 3 x 3 50 LUCRAREA 7 FILTRAREA IMAGINILOR 1 1 1 1 1 25 25 25 25 25 1 1 1 1 1 25 25 25 25 25 1 1 T 1 1 25 25 25 25 25 1 1 1 1 1 25 25 25 25 25 1 1 1 1 1 25 25 25 25 25 Figura 7 3: Mască de filtrare pătrată de dimensiune 5 x 5 1 1 1 1 1 1 1 49 49 49 49 49 49 49 1 1 1 1 1 1 1 49 49 49 49 49 49 49 1 1 1 1 1 1 1 49 49 49 49 49 49 49 1 1 1 1 1 1 1 49 49 49 49 49 49 49 1 1 1 1 1 1 1 49 49 49 49 49 49 49 1 1 1 1 1 1 1 49 49 49 49 49 49 49 1 1 1 1 1 1 1 49 49 49 49 49 49 49 Figura 7 4: Mască de filtrare patrată de dimensiune 7 x 7 (a) (b) Figura 7 5: Filtrarea de mediere: (a) imaginea originala; (b) imaginea filtrata cu o masca 7x7 Filtrul de mediere nu este folosit în practica, deoarece, pe lângă zgomot, are de suferit și semnalul util (vezi Figura 7 5) Filtrarea de mediere este deranjantă pentru imagine în acele zone în care imaginea conține frecvențe înalte (variații bruște), pentru că duce la apariția fenomenului de încețoșare1 1 (engl ) blurring 7 1 FILTRAREA LINIARA A IMAGINILOR 51 Din punctul de vedere al zgomotului, pentru filtrare este utilă o mască de filtrare de dimensiune cât mai mare Din punctul de vedere al semnalului util, al imaginii, este util ca masca sa fie cat mai mică în practică se realizează un compromis între cele două aspecte Alte măști de filtrare Alte măști de filtrare, pentru filterele de mediere, sunt prezentate în continuare, deși toate au același randament nesatisfăcător / -1 X A \ 16 16 16 2 4 2 16 16 16 121 16 16 16 1 16 1 16 116 0 1 5 0 0 1 8 0 1 1 x 1 1 1 1 1 1 16 2 16 5 5 5 8 4 8 1 16 1 16 116 0 1 5 0 J 0 1 8 0 7 1 2 Filtrele trece-sus Filtrele trece-sus urmăresc eliminarea componentelor de frecvenȘă joasă din imagine Sunt folosite în general pentru detectarea frontierelor sau contururilor din imagine, acolo unde au loc treceri sau variaȘii bruște ale luminanȘei Condișia de normare a coeficienȘilor filtrelor trece-sus CondiȘia de normare a coeficienȘilor pentru un filtru trece-jos (relaȘia (7 6)) se determină impunand condiȘia ca filtrul să rejecteze complet (sau să atenueze complet) componenta continuăa a imaginii Pentru aceasta vom considera, la fel, o imagine f (i,j) = g pentru Vi, j La ieȘirea filtrului trece-sus vom avea g(i,j) = 0 pentru Vi,j 0 Wkl' V (7 5) (k,i) ew XX Wkl = 0 (7 6) (k,i) ew Filtrul de accentuare Filtrul de accentuare nu este un filtru trece-sus, ci foloseȘste filtrarea trece-sus pentru a realiza accentuarea Prin accentuare se înȘtelege contrastarea unei imagini Și are ca scop îmbunătăȘirea perceperii vizuale a contururilor obiectelor Cu alte cuvinte, îmbunaătaăȘtirea detectabilitaăȘtii componentelor scenei de-a lungul frontierelor acestora Acest lucru se realizează, în principiu, prin modificarea valorilor pixelilor situaȘi de o parte Și de alta a unei frontiere comune 52 LUCRAREA 7 FILTRAREA IMAGINILOR Sistemul uman are tendința de a adânci profilul zonelor de tranziție dintre regiuni uniforme Studiul fiziologiei sistemului vizual a demonstrat că aceasta se realizează prin prelucrari de tip derivativ ce apar în diferitele etape pe care le parcurge informația vizuala Efectul global poate fi modelat matematic prin scăaderea din semnalul original a unei derivate secunde ponderate în continuare sunt prezentate cateva mățti de implementare a unei derivate secunde de tip Laplace: / 0 - 4 0 \ 1 i 1 4 1 4 \ 0 - 4 0 / / 1 1 1 \ ' 4 2 4 1i 1 2 1 2 l 1 -1 1 ) \ 4 2 4 / 1 1 1 ' 8 8 8 ' 1 i 1 8 8 1 1 1 \ 8 8 8 / în Figura 7 6 puteți observa efectele filtrarii Laplace (b) Figura 7 6: Filtrarea Laplace: (a) imaginea originală; (b) imaginea rezultată (negativată) Filtrul de accentuare se implementează după schema prezentată în Figura 7 7 7 2 Filtrarea neliniară a imaginilor Filtrele neliniare nu respectăa principiul liniarităațtii sau al superpozițtiei enunțat la începutul lucrarii Acestea au aparut din necesitatea de a depăsi limitările filtrelor liniare, în ceea ce privețte zgomotele care nu au o distribuție normalăa sau nu sunt aditive 7 2 FILTRAREA NELINIARA A IMAGINILOR 53 Figura 7 7: Filtrul de accentuare 7 2 1 Filtrele de ordine Filtrele de ordine sunt operatori locali, definiți la rândul lor de o fereastră, care selectează din imagine un numar de pixeli vecini pixelului curent, ântr-un mod identic cu tehnica ferestrei glisante Valorile pixelilor selectați se ordoneazăa crescaător Să presupunem că fereastra conține n pixeli, iar valorile lor formează urmaătoarea mulțtime: {x1,x2, ,xn} (7 7) Dupa ce aceste valori au fost ordonate crescator, vom avea: {X(1) ;X(2) , ,X(n) } (7 8) pentru care sunt îndeplinite urmatoarele condiții: x(1) — x(2) — — x(n) (7 9) Ieșirea filtrului de ordine de ordin k, pentru k 2 întreg, este statistica de ordinul k, cu alte cuvinte, elementul de pe poziția k din șirul ordonat: rankk{x1 ,x2 , ,xn} = X(k) (7 10) Filtrul median Filtrul median ete un filtru de ordine a căarui iețsire este statistica de ordin central, adică elementul ce se află pe poziția de mijloc a țirului ordonat de valori selectate de fereastra de filtrare: )X( n+i) daca n este impar, 1 k 2i , (7 11) 2X( n) + 2X( n +1) daca n este par 54 LUCRAREA 7 FILTRAREA IMAGINILOR Filtrul median este potrivit pentru eliminarea zgomotului de tip “sare și piper" După ordonarea valorilor pixelilor, valorile zgomotului (adică 0 sau 255) se vor situa pe primele, respectiv ultimele poziții în mulțime, si deci, la iețirea filtrului, vom avea o valoare diferită de valorile zgomotului Totuți există ți situații în care, după filtrare, mai exista pixeli afectați de zgomot In acest caz, spunem că filtrul a fost “strapuns" de zgomot Acest lucru este posibil atunci cand mai mult de jumătate din pixelii selectatți de fereastra de filtrare, sunt afectați în acelați mod (sare sau piper, 255 sau 0) de zgomot (a) (b) Figura 7 8: Filtrul median: (a) imaginea originală afectată de zgomot “sare si piper"; (b) imaginea filtrată Filtrul de minim Filtrul de minim este un filtru de ordine a căarui iețsire este statistica de ordinul 1, adică valoarea x(p, care este cea mai mica valoare din mulțimea de valori ale pixelilor selectațti de căatre fereastra de filtrare Filtrul de maxim Filtrul de maxim este un filtru de ordine a căarui iețsire este statistica de ordinul n, adică valoarea X(n), care este cea mai mare valoare din mulțimea de valori luate în considerate DESFĂȘURAREA LUCRĂRII Problema 1 Observați efectele filtrului de mediere pentru o imagine afectată de zgomot gaussian Codul C al funcției care implementează filtrul de 7 2 FILTRAREA NELINIARA A IMAGINILOR 55 mediere este prezentat în continuare: void ImageViewer :: filtru mediere( void ) { int i, j; int k, l; int w, h; double v[ 3 ][ 3 ]; //coeficientii mastii de filtrare v = 1 /9; v = 1 /9; v = 1 /9; v = 1 /9; v = 1 /9; v = 1 /9; v = 1 /9 v = 1 /9 v = 1 /9 w = image width(); h = image height(); QImage image fil( w, h, 32, 0, QImage::IgnoreEndian ); for( i = 1; i sir[ m + 1 ] ) { aux = sir[ m ]; sir[ m ] = sir[ m + 1 ] sir[ m + 1 ] = aux; k = 0; } 58 LUCRAREA 7 FILTRAREA IMAGINILOR } //elementul median med = sir[ 4 ]; //noua valoare a pixelului image fil setPixel( i, j, qRgb( med, med, med ) ); } image = image fil; pm = image; update(); } Problema 7 Implementați filtrul de minim Observați efectele lui asupra unei imagini neafectate de zgomot Problema 8 Implementați filtrul de maxim Observați efectele lui asupra unei imagini neafectate de zgomot Lucrarea 8 Transformări unitare BREVIAR TEORETIC Transformările reprezintă o categorie de prelucrări ce include operații de tip integral, la calculul noii valori a unui pixel al imaginii transformate contribuind valorile tuturor pixelilor din imaginea originală Termenul de transformare se referaă la o clasăa de matrici unitare folosite pentru a reprezenta imagini O matrice A de dimensiune N x N este unitară daca inversa ei este matricea A*T: A • A”1 = A • A*T = A*T • A = IN (8 1) unde * reprezinta operația de complementare în mulțimea numerelor complexe, T reprezintă operația de transpunere a unei matrici, iar In este matricea identitate de dimensiune N x N: 0 \ 0 0 1 0 0 0 1 0 0 0 1 In = \ 0 0 0 (8 2) 1 Pentru un vector uni-dimensional u de dimensiune N, de forma: u(0) u(1) u = = [u(0), u(1), , u(N — 1)]T (8 3) u(N - 1) se numețte transformare unitară directă relația: 59 60 LUCRAREA 8 TRANSFORMĂRI UNITARE N-1 v(k) = a(k, n) ■ u(n), 0 max ) max =C[i][j]; } //formarea matricei U for( i = 0; i max ) max =V[i][j]; } QImage transf( N, N, 32, 0, QImage::IgnoreEndian ); for( i = 0; i max ) max =V[i][j]; } //anularea coeficientilor //in vederea maririi factorului de compresie for( i = 0; i 1 1 1 ) 4 2 4 \ 1 1 1 / \ 8 8 8 / 1 8 1 1 8 Pentru fiecare pixel din imagine se va calcula derivata a doua ți dacă aceasta este zero, atunci pixelul este considerat ca fiind pixel de contur DESFĂȘURAREA LUCRĂRII Problema 1 Implementați segmentarea prin praguirea histogramei Pentru aceasta veti folosi codul de la lucrarea nr 4 (egalizarea histogramei) ți imaginea tools bmp Problema 2 Observați extragerea contururilor prin tehnici de gradient, folosind un operator Prewitt void ImageViewer :: extragere contururi( void ) { int i, j; int k, l; int w, h; 88 LUCRAREA 10 SEGMENTAREA IMAGINILOR int g, g1, g2; int h1[ 3 ][ 3 ], h2[ 3 ][ 3 ]; //operatori Prewitt h1 = -1; h1 = -1; h1 = -1; h1 = 0 h1 = 0 h1 = 0 h1 = 1; h1 = 1; h1 = 1; h2 = -1; h2 = -1; h2 = -1; h2 = 0; h2 = 0; h2 = 0; h2 = 1; h2 = 1; h2 = 1; w = image width(); h = image height(); QImage image ext( w, h, 32, 0, QImage::IgnoreEndian ); for( i = 1; i 100 ) image ext setPixel( i, j, qRgb(255,255,255) ); else image ext setPixel( i, j, qRgb(0,0,0) ); } image = image ext; pm = image; update(); } 89 10 2 SEGMENTAREA ORIENTATA PE CONTURURI Problema 3 Implementați operatorul Izotrop de extragere a contururilor Problema 4 Implementați operatorul Sobel de extragere a contururilor Problema 5 Implementați operatorul compas dat drept exemplu Problema 6 Implementați extragerea de contururi prin identificarea trecerilor prin zero ale derivatei secunde Problema 7 Implementați crețterea de regiuni, pornind de la un singur germene Veți folosi pentru aceasta imaginea tools bmp 90 LUCRAREA 10 SEGMENTAREA IMAGINILOR Lista figurilor 2 1 Formarea imaginii în aparatul de fotografiat: 1-lentilă, 2- camera obscura, 3-obiectiv, 4-peliculă 6 2 2 Forma tipică a funcției de eficiența luminoasă relativă 7 2 3 Domeniu finit din R2 8 2 4 Rețea dreptunghiulară de ețantionare 9 2 5 Exemplu de funție de cuantizare 10 3 1 Reprezentarea unei operații punctuale 16 3 2 Funcția de accentuare de contrast 16 3 3 Funcția de întindere maxima a contrastului 17 3 4 întinderea maximă a contrastului: (a) imaginea originală ți (b) imaginea rezultată, pentru a=50, b=200 18 3 5 Funcția de binarizare 19 3 6 Binarizarea: (a) imaginea originală ți (b) imaginea binarizată cu T=127 19 3 7 Funcția de negativare 20 3 8 Negativarea: (a) imaginea originala ți (b) negativul imaginii 20 3 9 Funcția de decupare cu păstrarea fundalului 21 3 10 Funcția de decupare fără pastrarea fundalului 21 4 1 Egalizarea histogramei: (a) imaginea originala ți (b) imaginea rezultatăa 28 4 2 Histograma (a) originală ți (b) dupa egalizare 28 5 1 Translatarea unui obiect dreptunghiular 32 5 2 Translatarea spre dreapta a unui obiect dreptunghiular într-o imagine 32 5 3 Rotația 33 5 4 Oglindirea (a) “stanga-dreapta” (b) “sus-jos” 33 5 5 Oglindirea “stanga-dreapta”: (a) imaginea originală; (b) imag- inea oglindită 34 6 1 Modelul aditiv de degradare 38 91 92 LISTA FIGURILOR 6 2 Funcția de densitate de probabilitate pentru o distribuție uniformă 39 6 3 Zgomotul uniform: (a) imaginea originală; (b) imaginea afectată de zgomot uniform, SNR=5 dB 39 6 4 Funcția de densitate de probabilitate pentru o distribuție gaussiană 40 6 5 Zgomotul Gaussian: (a) imaginea originala; (b) imaginea afectatăa de zgomot gaussian, SNR=5 dB 40 6 6 Zgomotul “sare ți piper": (a) imaginea originală; (b) imaginea cu 10% pixeli afectați de zgomot 41 6 7 Diverse funcții de distribuție 42 7 1 Mască de filtrare patrată de dimensiune 3x3 48 7 2 Masca de filtrare pătrata de dimensiune 3 x 3 49 7 3 Masca de filtrare patrată de dimensiune 5 x 5 50 7 4 Mascăa de filtrare paătrataă de dimensiune 7 x 7 50 7 5 Filtrarea de mediere: (a) imaginea originală; (b) imaginea filtratăa cu o mascăa 7x7 50 7 6 Filtrarea Laplace: (a) imaginea originală; (b) imaginea rezultată (negativată) 52 7 7 Filtrul de accentuare 53 7 8 Filtrul median: (a) imaginea originalăa afectatăa de zgomot “sare ți piper"; (b) imaginea filtrată 54 9 1 Algoritmul Huffman: etapa 1 71 9 2 Algoritmul Huffman: etapa 2 71 9 3 Algoritmul Huffman: etapa 3 72 9 4 Algoritmul Huffman: etapa 4 72 9 5 Transformarea imaginii într-un șir unidimensional, prin concatenarea liniilor 73 9 6 Transformarea unei imagini cu 256 nivele de gri, în 8 imagini binare 74 9 7 Transformarea directăa 75 9 8 Anularea coeficiențtilor de energie micăa 75 9 9 Transformarea inversăa 75 10 1 Exemplu teoretic de segmentare 82 10 2 Histogramă: (a) bimodală; (b) cu trei moduri dominante 83 10 3 Filtrul Sobel: (a) imaginea originală; (b) imaginea rezultată 86 Bibliografie V Buzuloiu, “Prelucrarea Imaginilor", note de curs C Vertan, “Prelucrarea și Analiza Imaginilor”, Ed Printech, 1999 A Jain, “Fundamentals of Digital Image Processing”, Prentice-Hall, 1989 R Gonzales, R Woods, “Digital Image Processing”, Addison Wes-ley, 1992 Trolltech, “Qt OverView”, http://www trolltech com/ Trolltech, “The Qt Class Reference”, http://www trolltech com/ 93 ISBN (10) 973-635-674-4 ISBN (13) 978-973-635-674-2